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Abstract. We study a critically coupled cavity doped with resonant atoms with 
metamaterial slabs as mirrors. We show how resonant atom-cavity interaction can 
lead to a splitting of the critical coupling dip. The results are explained in terms of 
the frequency and lifetime splitting of the coupled system. 
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1. Introduction 

Over the past decade very many exotic applications of negative index materials 
(NIMs), otherwise known as metamaterials or left-handed materials have been suggested 
and tested. These range from superlensing to lasing spasers, optical nano-circuits 
to invisibility cloaks, cavity quantum electrodynamics (QED) to electromagnetically 
induced transparency [H [21 [3l HI O [6] . In most of these applications the doubly negative 
response of the metamaterials (namely, negative permittivity e and permeability n) 
over a certain range of frequencies has been exploited. In fact, frequency dispersion 
and absorption have been identified as some of the crucial characteristics of the 
metamaterials, determining their figure of merit and suitability for device applications 
[TJ m [9l [TOl [11]. Only recently it has been recognized that metamaterials, capable 
of exhibiting negative index, can have interesting applications in a different frequency 
domain (where it is not truly a negative index material). It was shown that a slab 
of a metamaterial can exhibit a stop gap (as in one dimensional periodic structures) 
in the frequency range between the electric and magnetic plasma frequencies [T2] . 
Interesting features of such a material in the said frequency domain follows from close 
to zero real part of the refractive index. Depending on the separation between the 
plasma frequencies, one thus has a material with broadband absorption features. These 
properties were exploited to suggest a cavity with metamaterial layers as mirrors, which 
can be critically coupled to incident radiation, resonant with one or more of the modes 
of this cavity. Recall that a critically coupled system can absorb almost all the incident 
light, resulting in near-null reflection as well as transmission. Critical coupling (CC) in 
layered structures with absorbing polymer films or metal nano composite layers have 
been studied in detail [131 [HI [IS]- It is understood that there are certain advantages 
for critical coupling in a metamaterial cavity (like flexible tunability, etc.) as compared 
to the earlier structures. These advantages follow from the fact that now one uses the 
Fabry-Perot (FP) modes (tunable by changing the cavity free spectral range) without 
any constraints imposed by the narrowband feature of the polymer or the nano composite 
absorber. Moreover, there is a field enhancement in the cavity, whenever it is critically 
coupled. In fact, these are the special features one looks for to study the strong 
interaction regime of cavity QED. Recall that in the weak interaction regime, the photon 
emitted by the excited atom can escape the cavity, before it can interact back with the 
atom. In case of strong interaction, due to high finesse of the cavity, the photon has a 
larger life time, and there is a periodic exchange of energy between the atom and the 
cavity mode. It is thus interesting to study atom-cavity interaction in the context of a 
critically coupled cavity. In this paper, we study the manifestations of this interaction 
when the intracavity medium is doped with atoms resonant with the critical coupling 
frequency. We show that like in usual vacuum field Rabi splitting [16] , here also we have 
the splitting of the critical coupling dip into two, one corresponding to the symmetric 
mode and the other to the anti-symmetric mode. We solve the dispersion relation 
for obtaining the corresponding branches for both the real and the imaginary parts 
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Figure 1. Schematic of ttie layered structure. 



of the frequency. The real part gives the approximate location of the modes, while 
the imaginary part gives the corresponding decay rates. Finally, we study the field 
distributions corresponding to the excitation of these modes. 

The structure of the paper is as follows. In Section 2, we define our system and 
formulate the problem. Section 3 is devoted to numerical results pertaining to the 
frequency splitting in the total scattering data and also the roots of the dispersion 
equation. Finally, in Conclusions, we summarize the main results of this paper. 



2. Formulation 



Consider the system shown in figure [T] consisting of a spacer layer between two slabs of 
metamaterials. Let the permittivity and permeability of the negative index materials 
be given by a Lorentz type response as follows [1] 
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where fep (fmp) is the electric (magnetic) plasma frequency, feo {fmo) is the electric (mag- 
netic) resonance frequency and 7 is the decay rate (we assumed the same decay rate 
for both electric and magnetic resonances). A slab of such a metamaterial was shown 
to exhibit a stopgap (like in one dimensional photonic band gap structures) between its 
electric and magnetic plasma frequencies, characterized by mostly imaginary refractive 
index [12]. We will refer to such metamaterials still as NIMs, though we consider a 
frequency domain where it does not possess negative refraction. We use two slabs of 
such NIMs with distinct values of the plasma frequencies such that the stop band of one 
(NIMt, see figure [I]) lies well inside that of the other (NIM^). We assume the NIM slabs 
to be isotropic and to have a 3D structure (as in [12]) although the experiment of Shelby 
et al [17] used a 2D structure. It was shown recently that a dielectric cavity formed 
by two such NIM slabs exhibits critical coupling near the edges of the stop band of 
NIM^ [18]. The interplay of the FP resonances of the cavity with the stop band features 
of the NIM layers resulted in the critical coupling phenomenon. Further, varying the 
angle of incidence, which controls the width of the NIMs' stop gap, was used to achieve 
critical coupling at other frequencies. Recall that there is practically no reflection or 
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transmission from a critically coupled cavity resulting in an almost perfect absorption 
of the incident energy by the structure. In this study our goal is to investigate the ef- 
fects of cavity-atom interaction when the intracavity medium (i.e., the dielectric layer) 
is doped with resonant atoms. We assume the atomic medium to be nonmagnetic and 
its response to be given by the following dielectric function 

where eh is the dielectric constant of the host material, is proportional to the 
dopant density, /02 and 72 are the resonance frequency and the decay rate of the atom, 
respectively. Furthermore, the atomic resonance frequency is assumed to be close to 
the CC frequency of the undoped cavity. In the following section we demonstrate that 
such structures can lead to splitting of the CC dip for sufficient atom-cavity coupling 
determined by the dopant density. Such normal mode splittings have been studied 
in detail in the context of FP, modulated FP or spherical cavities (supporting the 
whispering gallery modes) EDI ED 122] . 

It is clear that the reflection and transmission profile for a system shown in figure [1] 
can be calculated using the standard characteristic matrix approach [231 1211 [25]. Both 
the refiection and transmission coefficients have a common denominator, the poles of 
which bear the information about the characteristic frequencies of the system. The 
corresponding dispersion relation can be written as [20] [25] 

D = (mil + mi2Pf)Pi + (^21 + m22Pf) = , (3) 

where (i, j = 1, 2) are the elements of the total characteristic matrix of the structure. 
Pi J — l^iJ ^iJ (f*^^ '^El polarization) and 9ij are the angles of incidence and 

emergence in the first and the last medium, respectively. The dispersion relation Q 
can be solved only for complex frequencies, which carries all the information about the 
split modes and their associated decay rates. In fact, the real part of the roots gives the 
locations of the split modes, while the imaginary part corresponds to the width of these 
resonances. It is understood that it is nontrivial to solve equation and obtain the 
different branches. We use a nonlinear root finder and present the results in the next 
section. 



3. Results 

In what follows we use scaled frequencies (in units of /o = lOGHz = c/Aq) and lengths 
(in units of Aq) to enable us to deal with dimensionless quantities (as in [12]) and also for 
reference. We use the same set of parameters as in the earlier work to achieve critical 
coupling with an undoped spacer layer [18]. The parameters for the lower NIM slab 
(NIMfe) are as follows: Up = 1-095, f^p = 1.28, Uo = 1-005, feo = 1-03, 7 = 0.001 and 
^3 = 5, while for NIM^ we take the same values, except for f^p = 1-14, f^p = 1.175 and 
di = 3. This choice of parameters ensures that the stop band of NIM^ lies inside that of 
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Figure 2. Splitting of the CC resonance at normal incidence for two dopant atom 
decay rates (a) 72 = 0.0008 and (b) 72 — 0.001 with densities fp =0.005 (dashed), 
0.01 (dash-dot). The result for the undoped spacer layer {fp — 0) is reproduced (solid 
curve) for reference [I^. Here di — 3, d2 ^ 7, = 5. /02 ~ 1.174 is near degenerate 
to the CC frequency (when fp — 0). Note that all frequencies are scaled in units of /q. 



NIMfo. The host material for the spacer layer has e/j = 3.8. The entire layered structure is 
embedded in air (e^ = e/ = 1, /Xj = /i/ = 1). The spacer layer thickness (0^2 ~ 7) is fixed 
to attain CC for normal incidence near the right edge of the stop band of NIMj (solid 
curves in figure [2]) without doping the cavity {fp = 0) [18j. The resonance frequency of 
the dopant atom is chosen to be near-degenerate (not identical) with this CC frequency 
at /02 ~ 1.1738 to ensure a strong interaction of the atom with the cavity resonance 
at CC. We study the total scattering {R + T) of the structure at normal incidence 
using the dopant density as a parameter. The calculations for total scattering have 
been performed for two dopant densities, namely, fp = 0.005,0.01. The results for two 
decay rates (72 = 0.0008,0.001) are shown in figures Et^a) and (b), respectively. These 
figures clearly demonstrate the splitting of the CC dip as a manifestation of the strong 
atom-cavity interaction (compare with the results of references [201 [22] ). Unfortunately 
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Figure 3. CC is achieved (solid curves) for the split resonance at higher frequency (/) 
for two dopant decay rates (a) 72 = 0.0008, fp = 0.0265, / = 1.179 and (b) 72 = 0.001, 
fp = 0.0295, / = 1.18. The split resonances at fp — 0.01 (dashed curves) are produced 
for reference. The other parameters are the same as in figure [21 
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the critical coupling is lost at both the split frequencies. However, the higher frequency 
dip is closer to critical coupling for larger densities. Thus a larger dopant density can 
restore critical coupling at the right frequency component for a denser atomic system. 
This restoration is shown in figures [3]^a) and (b) for both the values of 72. For example, 
for 72 = 0.0008 (72 = 0.001), CC is restored at fp = 0.0265 (/p = 0.0295). The dopant 
density may thus be used as a parameter to tune the CC frequency of the coupled 
system. A comparison of the left and right panels of figures 2 and 3 reveal that the 
splittings are more in case of lower atomic decay, which is quite expected. In fact the 
resolution of the split resonances depend on how the coupling compares with the decay 
rates. The lower the decay rates, the higher should be the splitting. Another pertinent 
feature that should be noted from these figures is that the width of the higher (lower) 
frequency split resonance increases (decreases) with increasing fp. 

We now show that all the above features can easily be understood based on an 
analysis of the roots of the dispersion relation (|3]). The real and imaginary parts of the 
roots (as functions of the dopant density) for the two values of 72 (72 = 0.0008, 0.001) 
are shown in the left and right panels of figure H] . The solid and dashed lines give the 
two branches for both. The well known frequency and lifetime splittings [201 122] can 
easily be read from the top and bottom panels of figure |H It is clear that the split mode 
frequencies move away from each other with increasing fp (top panel of figure H]) in 
conformity with the previous total scattering results. The split resonance at the higher 
frequency (dashed curves) widens, whereas, the one at lower frequency (solid curves) 
becomes narrower with increasing fp (see the bottom panels in figure H]). 

The correspondence of the split resonances to symmetric and anti-symmetric modes, 
as in vacuum Rabi splitting, can be inferred from a study of the field distribution inside 
the structure. The absolute value of the field amplitude must reach zero inside the cavity 
for the antisymmetric mode but not for the symmetric one. We use the parameters 
in figure Et^a) and plot in figure |5] the absolute value of the amplitude at both the left 
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Figure 4. The real ((a) and (c)) and imaginary ((b) and (d)) parts of the roots of the 
dispersion relation (/r) of the layered structure as a function of the atomic density at 
7 — 0.0008 (left panel) and 0.001 (right panel). The solid (dashed) curve corresponds 
to the split resonance at the lower (higher) frequency. 




Signature of strong atom-cavity interaction on critical coupling 



7 




Figure 5. The field distribution inside the structure for the parameters in figure [3Ia). 
(a) fp — 0.01 with split resonances at / = 1.172 and 1.176 (b) fp = 0.0265 with 
split resonances at f — 1.168 and 1.179 . The solid(dashed) curves correspond to the 
lower(higher) frequency split resonance, z — corresponds to the air-NIMf interface 
and z = di + 1^2 + da ~ 15 to the NIMf,-air interface. The insets are enlarged portions 
of the corresponding plots. 

(antisymmetric) and right (symmetric) resonances (solid and dashed hnes, respectively). 
z in figure [5] corresponds to the depth inside the layered structure from the air-NIM^ 
interface. From the insets in figure |5] it is seen that the magnitude of the field amplitude 
touches zero for the antisymmetric mode only. Further, we note the significant field 
enhancement in case of the critically coupled cavity (see figure [5b) 

4. Conclusion 

In conclusion, we have studied a critically coupled cavity and shown that the critical 
coupling dip can exhibit a normal mode splitting when the intracavity medium is 
doped with resonant atoms. The intracavity medium with dopant atoms was modelled 
by a Lorentz type dielectric with characteristic resonance frequency and decay rate. 
Numerical results for the total scattering R+T was obtained using a characteristic 
matrix formulation. The corresponding dispersion relations were solved for the location 
and width of the split resonances. Moreover, it was demonstrated that the atomic 
dopant density can be used as an additional handle to achieve critical coupling with 
respect to one of the split resonances. 

References 

[1] Shalaev V M 2007 Nat. Phot. 1 41 

[2] Pendry J B 2000 Phys. Rev. Lett. 85 3966 

[3] Zheludev N I, Prosvirnin S L, Papasimakis N and Fedotov V A 2008 Nat. Phot. 2 351 
[4] Engheta N 2007 Science 317 1698 

[5] Xu J P, Yang Y P, Lin Q and Zhu S Y 2009 Phys. Rev. A 79 043812 

[6] Papasimakis N, Fedotov V A, Zheludev N I 2008 Phys. Rev. Lett. 101 253903 

[7] Veselago V G 1968 Soviet Physics Uspekhi 10 509 



Signature of strong atom-cavity interaction on critical coupling 



8 



[8] Pendry J B 2004 Contemporary Physics 45 191 

[9] Kinsler P and McCall M W 2008 Phys. Rev. Lett. 101 167401 

[10] Nistad B and Skaar J 2008 Phys. Rev. E 78 036603 

[11] Golla D, Deb S, Dutta Gupta S 2009 Competition between structural and intrinsic dispersion in 

delay through a left handed medium Preprint arXiv:0911.5374v2 [physics. optics] 
[12] Bloemer M, D'Aguanno G and Scalora M 2005 App. Phys. Lett. 87 261921 
[13] Tischler J R, Bradley M S and Bulovic V 2006 Opt. Lett. 31 2045 
[14] Dutta Gupta S 2007 Opt. Lett. 32 1483 

[15] Deb S, Dutta Gupta S, Banerji J and Dutta Gupta S 2007 J. Opt. A: Pure Appl. Opt. 9 555 
[16] Carmichael H J, Tian L and Ren W Advances in Atomic, Molecular and Optical Physics supplement 

2 ed P R Berman (Academic Press, 1994) 
[17] Shelby R A, Smith D R and Schultz S 2001 Science 292 77 

[18] Deb S and Dutta Gupta S 2009 Critical coupling in a system of two coupled negative index medium 

layers Preprint arXiv:0912.0575vl [physics. optics] 
[19] Manga Rao V S C and Dutta Gupta S and Agarwal G S 2004 Opt. Lett. 29 307 
[20] Dutta Gupta S and Agarwal G S 1993 Opt. Comm. 103 122 
[21] Agarwal G S and Dutta Gupta S 1992 Opt. Comm. 93 173 
[22] Dutta Gupta S and Agarwal G S 1995 Opt. Comm. 115 597 
[23] Yeh P Optical Waves in Layered Media (Wiley, 1988) 

[24] Born M and Wolf E Principles of Optics: Electromagnetic Theory of Propagation (Cambridge 
University Press, 2001) 

[25] Dutta Gupta S Progress in Optics {Nonlinear optics of stratified media vol 38) ed Wolf E (Elsevier 
Science, 1998) 



